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Sudden changes in the dynamical properties of a supercritical fluid model have been found as a function of 
pressure and temperature (T/Tj — 2-5 and P/P^ — 10-10'), striving with the notion of a single phase beyond 
the critical point established by thermodynamics. The sound propagation in the Terahertz frequency region 
reveals a sharp dynamic crossover between the gas like and the hquid like regimes along several isotherms, 
which involves, at sufficienfly low densities, the interplay between purely acoustic waves and heat waves. 
Such a crossover allows one to determine a dynamic line in the phase diagram which exhibits a very tight 
correlation with a number of thermodynamic observables, showing that the supercritical state is remarkably 
more complex than thought so far. 

The supercritical fluid has since ever', until recently, been considered as a unique state with intermediate 
properties of the liquid and gas phases. The thermodynamical observables show a continuous evolution 
along any P-T path in supercritical conditions. This fact is supported by the traditional view on liquid and 
gas as phases of the same symmetry, implying that these phases cannot be distinguished above the critical point 
where they do not coexist. On the other hand, this view assumes that physics of fluids is entirely determined by 
their structure. Modern approaches to liquid theory clearly show that this is not always the case. Dynamics plays 
often a crucial role. Early in 2005 L. Xu and collaborators have pointed out that the diffusion coefficient of water 
undergoes a non-Arrhenius ("fragile") to Arrhenius ("strong") crossover accompanying a thermodynamic 
liquid-liquid transition in the so called "no-man's land" i.e. not accessible by current experimental techniques^. 

More recently, it has been discovered that subtle but remarkable changes in supercritical fluids do occur in the 
positive sound dispersion (PSD)'"^. These changes posses the potential to part the phase diagram in a gas like and 
a liquid like region. 

PSD is a phenomenon ascribed to the presence of one or more relaxation processes consisting in the increase of 
the sound velocity when the sound frequency is close to the relaxation frequency. In fact the relaxation processes 
time scale, x, marks the border between the "viscous liquid-like" dynamic (for (nL(Q)T <C 1, where (Ol and Q are the 
sound frequency and wave vector, respectively) to the "elastic solid-like" one (col(Q)t3>1). Correspondingly, the 
longitudinal sound velocity undergoes a transition (PSD) from its "low" frequency adiabatic limit Cs to its 
"infinite" frequency limit c„ > Cg. This effect has been observed in a wealth of liquids and glasses with typical 
relaxation frequency of the order of THz in simple liquids. On the other hand, PSD has never been observed so far 
in gases'"'. Recently, PSD has been observed in supercritical fluids such as Argon (critical parameters T^ = 
150.87 K and P^ = 4.898 MPa) and Oxygen (T^ = 154.6 K and Pc = 5.050 MPa) along isothermal lines, showing 
a cusp (discontinuity in the first derivative) as a function of pressure. This allows to define the high and low 
density supercritical fluids as liquid-like and gas like supercritical fluids, respectively' '. Based on such results, it 
has been proposed that the onset of PSD takes place in correspondence with the maxima of constant pressure 
specific heat Cp, defining the so called Widom line', which establishes a tight link between dynamics (PSD) and 
thermodynamics' (cp). Unfortunately, in that case the study was performed along an isotherm at high temper- 
ature (T/Tj = 3.8), where the maximum of Cp is no longer visible and the Widom line could be only extrapolated. 
A recent study indeed identified a dynamic crossover at the Widom line at lower temperatures where the maxima 
in Cp are still very well visible, based on anomalies in the self-diffusion coefficient and the excess two body 
entropy". These results are at ease with previous findings of our group on the dynamic crossover related to 
collective modes. On the other hand, our view has been challenged by recent theoretical investigations. In one of 
these works' it has been found that the extrema of different thermodynamic functions rapidly smear out above the 
critical point; this led the authors to question the potential link between dynamics and thermodynamics in the 
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supercritical state, and to suggest the separation between liquid-like 
and gas-like regions to be rather determined by dynamics alone. The 
crossover in PSD of a model supercritical fluid was indeed assigned to 
a universal dynamic (Frenkel) line, where the equality between the 
relaxation time and the minimum period of transverse quasiharmo- 
nic waves holds'". Unfortunately, PSD was not calculated, hence the 
advanced relationship between the potential dynamic crossover and 
the Frenkel line remains at the level of an interesting speculation 
calling for further investigations. Another possible way to identify 
two regions in the supercritical phase builds on the appearance of 
some elasticity upon crossing an aggregation line". In addition, an 
alternative definition of the Widom line in terms of the loci of max- 
ima of the volume of the correlation length was recently advanced'^. 
Interestingly, this line does not follow the slope of any response 
function extrema except in the close vicinity of the critical point, 
where it indeed follows the Cp maxima. This result adds a further call 
for a definitive clue on what is the right thermodynamic observable 
-if any- to correlate with the dynamic crossover. 

Results 

AU this matter motivated us to deeply investigate what we think to be a 
milestone in the physics of the elusive supercritical state, i.e. the poten- 
tial link between dynamics and thermodynamics. Specifically, based on 
molecular dynamics (MD) computer simulations, we investigated the 
onset of PSD along several isotherms down to temperatures where the 
line of maxima of Cp is well defined. In fact, we discovered that the onset 
of PSD occurs along this line up to temperatures where the maximum 
of Cp is still well defined and along the common minima of thermal 
difiusivity, Dx, and kinematic shear viscosity, v, at higher temperatures 
(note that kinematic longitudinal and shear viscosity Dl and v are 
related to the shear and bulk viscosity coefficients r|s and rjg as: Dl 
= (4/3ris -I- r|B)/p> and v = ris/p> where p is the density). The lines of 
extrema of these three observables are very close each others at suffi- 
ciently low temperatures, where maxima of Cp are well defined. 
Remarkably, we found that at sufficiently low densities in the gas-like 
regime, up to the dynamic crossover to the liquid-like regime, there is a 
strong interplay between purely acoustic waves and heat waves. These 
findings led us to the quantitative characterization of both types of wave 
properties, along with the quantitative determination of the crossover. 

Some S(Q, (o)s simulated on Ar at 700 K at four selected pressure 
and Q values are reported in figure 1, upper panel. The acoustic peak 
is increasingly more pronounced at high densities and small Q 
values. The corresponding sound dispersion Cl(Q) values have been 
divided by the adiabatic sound velocity value Cs (figure 1, lower 
panel). In case of no PSD this ratio must be equal to 1 at small Qs 
and to monotonically decrease at higher Qs on approaching the Q 
value of the first sharp diffraction peak (FSDP) of the static structure 
factor S(Q), resembling the bending of dispersion curves in crystals 
on approaching the BriUouin zone boundaries. On the other hand, 
when PSD is present, as it is the case here, the Cl(Q)/cs ratio is 
expected to tend to 1 at very small Qs, to show a maximum upon 
increasing Q and then again to decrease at higher Qs on approaching 
the FSDP. PSD, i.e. the value of this maximum has been estimated by 
fitting a cubic polynomial function to the Cl(Q)/cs ratio to account 
for finite Q sampling. The maximum takes place at around 
3-5 nm~', i.e. at the lengthscale ~2n/Q*, where the inverse of struc- 
tural relaxation time t(Q) equals the acoustic frequency fi(Q). This 
occurs at all the investigated pressures and the amount of PSD 
decreases on lowering density. Similar evolutions have been obtained 
along the other isotherms. 

Density evolution of PSD is reported in Figure 2. Along all iso- 
therms PSD shows a kink, being constant at low densities and 
increasing linearly at high densities. The extrapolation of the linear 
behaviour to low densities allows to estimate the density at which 
PSD is expected to vanish if coupling to the structural relaxation 
process was the only underlying, relevant mechanism; this density 



marks the onset of this kind of PSD or, in other words, the dynamic 
crossover between the gas-like and the liquid-like regime. On the 
other hand, the non-vanishing plateau in PSD at low densities indi- 
cates that for p^O some mechanism other than structural relaxation 
must exist. In order to shed light on potential mechanisms of PSD we 
used a thermo-viscoelastic model considering the possible propaga- 
tion of heat and sound waves and the coupling between them". A 
powerful method in the framework of generalized hydrodynamics is 
the approach of Generalized Collective Modes (GCM), representing 
acoustic excitations coupled to non-hydrodynamic processes as 
dynamic eigenmodes" '^. The GCM approach explained correctly 
the "fast sound" phenomenon in disparate mass binary liquid alloys 
and predicted conditions of existence of non-hydrodynamic excita- 
tions like optic-like modes, heat waves etc""'". In Figure 3 (left 
frames) we show the frequency of acoustic modes and heat waves, 
as functions of Qs, at four selected densities at T = 363 K. Heat waves 
were known to exist in solids as "second sound" excitations"*"^" whOe 
in liquids heat waves similarly to transverse sound waves cannot 
propagate on macroscopic distances^'. This fact is reflected in the 
Q dependence of the frequency of heat waves vanishing at some finite 
Qh value; heat waves can exist on spatial scales not larger than 
roughly 27t/Qi,. The dependence of sound dispersion on density 
shown in Fig. 3 gives evidence that below 5 nm"' the deviation from 
the hydrodynamic dispersion law reduces with decreasing density. 
For low densities, however, there is also a strong increase of the 
deviation from the hydrodynamic dispersion around 10 nm"', that 
is responsible of the overall PSD causing the plateau observed in 
Figure 2. We suggest that the origin of this increase at around 
10 nm"' for low-density states may be connected with heat waves 
because in this Q region at low density we observed the strong 
increase of their contributions to the S(Q,co) (Figure 3, right frames). 
These contributions (mode strengths) were obtained fit-free from the 
eigenvectors associated with acoustic and heat excitations. In the 
long-wavelength limit they tend to the values predicted by hydro- 
dynamics^^. At high density the strength of the contribution of heat 
waves to the S(Q,co) is negligible, whereas it strongly increases upon 
decreasing density. This implies that at low densities, beyond the 
hydrodynamic region there exist two branches of strongly coupled 
acoustic -heat excitations, which in turn causes strong deviation from 
the hydrodynamic dispersion for the high-frequency branch. 

Discussion 

The reduced phase diagram (P/P^-T/Tc) for Argon is reported in 
figure 4. It shows the simulation points (open dots) and the points 
where structural relaxation related PSD vanishes (full stars), accord- 
ing to the linear extrapolation of high density data as discussed 
above. Lines emanating from the critical point have been obtained 
from NIST data^': the red, blue and grey lines are obtained from 
extrema of Cp, v and Dj, respectively. These are among all possible 
lines emanating from the critical point and which can be considered 
as extensions of the liquid-vapour coexistence line into the super- 
critical fluid region. The line of maxima of Cp can be traced up to 
about 450 K (T/T^ = 2.98), as above this temperature these maxima 
severely smear out. On the other hand, the minima of Dj and v are 
still observable, even if broad, up to the highest temperature of this 
study (700 K, T/T^ = 4.64). The lines of extrema of Cp, D-j- and v are 
very close each other. Stars indicating the onset of PSD lie very close 
to these lines, putting the link between the dynamic crossover and 
thermodynamics on a very solid background. Indeed, the very exist- 
ence of such a relationship has been the matter of what now appears 
to be a correct intuition in our previous works'*'^, where the relevant 
marker was thought to be only the line of Cp maxima. On the other 
hand, the fact that this line is not defined far fi^om the critical point, 
where Cp maxima are not visible, made the connection between the 
onset of PSD and thermodynamics rather limited. Now, the connection 
between the dynamic crossover and the lines of D-p and v minima, 
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Figure 1 | Upper panel: S(Q, co) spectra obtained from MD simulations on supercritical Ar at four selected pressure (density) and Q values, at 700 K. 
Lower panel: Cl(Q)/cs ratio at the same pressures. The dispersions (dots) have been fitted to a cubic polynom for a better estimation of the maximum. 
The Cs values have been obtained from NIST23 up to the density corresponding to 10 Kbar, and extrapolated above this density using the formula: 
Cs(m/s) = 501.1174 + 0.1 7806 •p(kg/m") + 4.32538(10-''-p(kg/m")l We propose this extrapolation as, on the basis of NIST" data for Ar at 700 K, it 
results that in the whole range of available data (from 50 to 10000 bar, (33.866 to 1506.5 kg/m') the sound speed shows a quadratic dependence on density 
with very good approximation (within about 0.1%). 



which are well defined also very far from the critical point, makes the 
link between dynamics and thermodynamics quite robust and general. 
Interestingly, the non monotonic density dependence of longitudinal 
kinematic viscosity, Dj and v has a remarkable consequence. These 
quantities define the hydrodynamic relaxation times for longitudinal 
and transverse dynamics of supercritical fluids^^. The Q-width of the 
hydrodynamic regime, where collective excitations and relaxation pro- 
cesses have strictiy hydrodynamic behaviour, depends on inverse of 
kinematic viscosities and D-y^*-^^. Recentiy the estimates for the hydro- 
dynamic Q-range were proved by IXS experiments^*". Hence, the min- 



imum in the density dependence of these quantities corresponds to the 
widest Q-range of the hydrodynamic regime at given temperature, 
which in turn we foimd corresponding to the dynamic crossover 
observed in this work. Within this context, one may wonder how 
the crossover in the positive sound dispersion, established by this work 
to be the key parameter which allows to identify two different "phases", 
one gas- and the other liquid-like, relates to the behaviour of diffusion 
coefficient, playing a role in distinguishing two different liquid phases 
in water'. In the water case, Kumar et al^^'^". identify the Uquid-Uquid 
transition reporting the breakdown of a linear dependence of the dif- 
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Figure 2 | Density evolution (dots) of the PSD (max Cl(Q)/cs) at the four 
isotherms of this study. A Hnear behaviour has been fitted through the high 
density data points, whose intercept with the horizontal line of value 1 marks 
the onset of PSD, i.e. the dynamic crossover between the gas-like and the 
liquid-like regimes. The other Unes have been obtained using NIST data for 
Argorf' (corresponding y axes are arbitrarily rescaled) and represent: Cp 
(red), Dx (grey) and v (blue). The extrema of these Hnes correspond to the 
dynamic crossover marked by the onset of PSD. 

fiisivity as a function of the inverse structural relaxation time l/Xg,. By 
virtue of the MaxweU relation, in turn, they eventually put forward the 
breakdown of the Stokes-Einstein relation. Although fascinating, this 
latter step relies on the assumption of a constant "infinite shear modu- 
lus", the proportionality factor between diffusivity and l/Xo,. The pre- 
sent work addresses the non constant behaviour of positive dispersion 
as function of the thermodynamic coordinates in the case of the gas- 
liquid transition in supercritical fluids, precisely the deviation of infinite 
frequency elastic properties fi-om their hydrodynamic limit. It strongly 
suggests, therefore, that the amplitude of the structural relaxation was 
the dynamical observable related to the Widom line, rather than the 
relaxation time itself This observation seems to indicate that the two 
phenomenon (Uquid-liquid-transition in water and gas-Uquid trans- 
ition in supercritical fluids), both related to the Widom line, may be not 
based on the same physics. In summary, we discovered two striking 
features in density dependence of PSD in supercritical fluids: i) linear 
behaviour of PSD at sufficiently high densities, that extrapolates right 
to the density of Cp, D-p and v extrema; ii) a plateau in the density 
behaviour of PSD for clearly indicating that this observable is 

originated by an additional collective process along with structural 
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Figure 3 | Dispersion and mode strengths in S(Q, (b)s of two types of 
collective excitations in supercritical Ar at T = 363 K as follows from the 
fit-free analysis based on GCM approach. Left frames: dispersion of 
acoustic excitations (red) and non-hydrodynamic heat waves (green) at 
four selected densities. Dashed lines correspond to the hydrodynamic 
dispersion. Right frames: mode strengths of the corresponding 
propagating modes. In the long-wavelength limit heat waves do not exist, 
while the hydrodynamic value for the sound mode strength is 1/y, where 
the ratio of specific heats y varies from 1.43 to 2.16 for the highest and the 
lowest density, respectively. 

relaxation. The detailed analysis of dynamic eigenmodes and their 
contributions to S(Q,a))s revealed a novel effect of non-hydrodynamic 
heat waves on sound dispersion in low-density supercritical fluids. We 
quantitatively identified the dynamic crossover in a supercritical fluid 
model, at several temperatures, between the low density gas-like regime 
and the high density Uquid-hke regime, based on the analysis of col- 
lective modes with mesoscopic/microscopic wavelengths. In this 
framework, the unique role of heat thermal waves along with their 
interplay with sound waves have been unveiled. A remarkable outcome 
of this study is that the dynamic crossover shows a very tight correla- 
tion with the extrema of thermodynamic response functions and trans- 
port coefficients as Cp, Dt and v, defining a narrow region emanating 
from the critical point, as emphasized in reference^'-". This gives a new 
dimension to a recent work where a crossover in single particle 
dynamics was identified at the line of maxima of Cp in a supercritical 
model fluid", in turns anticipated by similar results in the different 
context of the liquid-liquid transition in water. The deep link between 
dynamics and thermodynamics beyond a critical point is now ratio- 
nalized on a very general background. 

Methods 

Molecular dynamics simulations were carried out on supercritical Argon at 
different densities at T — 363, 453, 573 and 700 K. The corresponding pressures were 
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Figure 4 | Reduced phase diagram (P/Pc -T/Tc) for Argon. Open dots: P-T 
points used in MD simulations. Stars: onset of PSD at the four isotherms, 
based on Hnear extrapolation of the PSD data points to low densities (see 
figure 2). Lines (obtained from NIST"): loci of the extrema of Cp (red), D^ 
(grey) and v (blue). Dotted line: critical isochore. 

obtained from NIST^^ up to 1 GPa and, above this pressure, as: n(g/cm^) = A + 
B-log(P(GPa)), where A - -1.41806 g/cm^ and B - 9.42809-10-* g/cm^ at 700 K, 
while the values at 573 K are reported elsewhere^. MD was done on the standard 
microcanonical ensemble for 2000 particles interacting through an ab initio poten- 
j.^gpo,3i Y2 A cutoff radius. We checked, at a few densities, that the same results are 
obtained with 4000 particles and 20 A cutoff radius. The production runs were of 
600000 time steps, for small densities they were extended up to 1.2 mln time steps, 
each one of 2 fs. The energy drift over each production run was not higher than 0.2% 
that gave evidence of very good energy conservation in MD simulations. All of the 
calculated thermodynamic quantities agree with the values provided by the NIST^^ 
within 1%. Transport coefficients reported in Fig. 2, such as D p and v were taken from 
the same NIST^' database. 

Every sixth timestep we sampled spatial Fourier transforms of density of particles 
n(Q,t), density of mass-current J(Q,t), and energy density e(Q,t) from the instantaneous 
positions ri(t), velocities Vi(t) and single -particle energies Si(t) of N particles as follows: 

These densities represent fluctuations of conserved quantities and are the main 
dynamic variables for hydrodynamic description of liquid dynamics. The density- 
density time correlation functions F„n(Q,f) — < n* {QJS)n{Q,t) > , where the asterisk 
means complex conjugation, were calculated from the MD time evolution of n(Q,t), 
and additional averages were performed over the all possible directions of wave 
vectors corresponding to the same wave number Q. Numerical Fourier transform was 
used to obtain from the calculated density- density time correlation functions the 
dynamic structure factors S(Q,co) and longitudinal current spectral functions Jl(Q,<b) 
= ((o/Q)^S(Q,co). The longitudinal sound velocity Cl(Q) is obtained as Cl(Q) ~ 
roL(Q)/Q) where (i)l(Q) is the sound frequency and it corresponds to the frequency of 
the maximum of Jl(Q)W)^- For analysis of contributions from collective modes to 
dynamic structure factors we applied the approach of Generalized Collective Modes 
(GCM)^*''^ , that consists in solving the matrix form of generalized Langevin equation 
in terms of dynamic eigenmodes that can exist in the system. The MD-derived time 
evolution of dynamic variables was used to calculate for each wave number the matrix 
elements of generalized hydrodynamic matrix T(Q)^*'^^. The eigenvectors associated 
with each eigenvalue were used to calculate the Q-dependent contributions (mode 
strengths) to dynamic structure factors. 
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